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Abstract 

In this paper we establish different representations of the so-called Yor integral, 
which is one of the key ingredient in mathematical finance, in particular, to compute 
normalized prices of Asian options. We show, that the Yor integral is related with 
the Kontorovich-Lebedev transform. Also we discuss its relationship with a system 
of polynomials recently introduced by the author. We derive new important prop- 
erties of these polynomials, including upper bounds, an exact asymptotic behavior 
for large values of their degree and explicit formula of coefficients. 
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1 Introduction and preliminary results 

In 1980 Yor [14] (see also in [3], [4], [15]) expressed the density of the Hartman -Watson 
distribution [2], which is related to the pricing of Asian options in mathematical finance, 
in the form of elementary integral, involving exponential functions 

e n2/2t f°° ( y 2 \ /7ry\ 
F t (r) — —== / exp I I exp (— r coshy) sinhysin ( — j dy, r, £ > 0. (1-1) 

However, despite of the importance of integral (1.1) in applications and a necessity to 
calculate it in the closed form, this task is quite difficult. Nevertheless, we can perhaps 
stimulate the corresponding numerical calculations of the integral and finding its asymp- 
totic behavior for small values of t, representing (1.1) in a different form. Precisely, we 
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will show the relation between F t (r) and the Kontorovich-Lebedev transformation [9], 
[10], [11] 

/•oo 

(G/)(r)= / K iT (r)f(r)dr, r e M+, (1.2) 
./o 

where K iT {r) is the modified Bessel function of the pure imaginary index ir [1]. As it is 
known, operator (1.2) is bounded 

G : L 2 (R+; rrfr) ^ L 2 (1R + ; r sinh vrrrfr) 

and the integral in (1.2) converges with respect to the norm in L 2 (R+; r sinh nrdr). More- 
over, it forms an isometric isomorphism between these Hilbert spaces and the Parseval 
identity holds 

/oo 2 roc 

Tsmhnr\(Gf)(T)\ 2 dT = Y j^ \f(r)\ 2 rdr. (1.3) 

Reciprocally, the inversion formula 

2 r°° 

f(r) = — / rsmhnrK tT (r)(Gf)(r)dr, (1.4) 

takes place, where the convergence of the integral (1.4) is understood with respect to the 
norm of the space L 2 rdr). 

The modified Bessel function K ir (x) is an eigenfunction of the following second order 
differential operator 

A x = x 2 -x— x— , (1.5) 

i.e. we have 

Ax K lT {x) = r 2 K iT {x). (1.6) 
It has the asymptotic behavior (cf. [1] relations (9.6.8), (9.6.9), (9.7.2)) 

1 /2 

= e- 2 [l + 0(l/^)], z^oo, (1.7) 

and near the origin 

#„(*) = O (z^ Rcu \) , z ->■ 0, (1.8) 

i^o(z) = -log* + 0(l), 2^0. (1.9) 

Moreover it can be defined by the following integral representations [9, (6-1-2)], [8], Vol. 
I, relation (2.4.18.4) 

PCX) 

K u (x)= e- xcoshu coshvudu,x>0, (1.10) 
Jo 



'0 

1 (X\ v 



cc 



,2 



K v {x) = 2 (i ) / e'^t-^dt, x>0. (1.11) 



o 
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The convolution operator for the Kontorovich-Lebedev transform is defined as follows [10, 
11] 



i poo poo 1 I x I yu I 

(f*h)(x) = (f(x)*g(x)) = —J J e 2 V - * J f(u)h(y)dudy, x > 0. (1.12) 

It is well defined in the Banach ring L a (R + ) = Li(R + ; K a (x)dx), «GR, i.e. the space of 
all summable functions / : R+ — > C with respect to the measure K a (x)dx for which 

POO 

||/||l«(r +) = / \f(x)\K a (x)dx (1.13) 
is finite. The following embeddings take place 



L Q (R+) = L~ a (R+), L a (R + ) C i7*(R+), |a| > > 0,a,/3 G R, 

2 

L a (R) D L P (R + ; xda:), 2 < p < oo, \a\ < 1 - -, 
where L p (R + ; xcix) is a weighted Banach space with the norm 

/ POO \ 1/p 

l p (r + ;^) = f y |/(x)| p a;rfxj , 1 < p < oo, 

H/IUooCR+jxdr) = ess sup x6R+ |/(a:)|. 

The factorization property is true for the convolution (1.12) in terms of the Kontorovich- 
Lebedev transform (1.2) in the space L Q (R + ), namely 

(G[f * h])(r) = (Gf)(r)(Gh)(r), r G R + . (1.14) 

This property is based on the Macdonald formula [1] 



1 r°° 1 U x +y i xy ) dt 
K v (x)KM = \\ q e 2 V - +t )K v (t)j. (1.15) 

2 A system of polynomials 

In this section we will provide a useful information about a system of polynomials, which 
is related to the Kontorovich-Lebedev transformation (1.2) and was studied by the author 
in [12]. In particular, we will derive new properties of these polynomials, including an 
upper bound, a series representation by the index of the modified Bessel functions of the 
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third kind and an explicit formula expressing their coefficients. In fact, as it is proved in 
[12], the following functions 

p n (x) = (-l) n e x AZe- x , neN , (2.1) 

where A™ is the n-th iteration of the differential operator (1.5), are n-th degree polyno- 
mials, which have the following integral representation 

2(—l) n f°° 

p n (x) = - -e x / r 2n K iT (x)dr, x > 0, n G N. (2.2) 

* Jo 

Hence calling the reciprocal formula (1.2), we obtain the Kontorovich-Lebedev transform 
of p n (x) 

/ K iT (x)e-* Pn (x)- = (-1)" — rr- y (2.3) 
J x smh(7rr) 

The system p n satisfies the differential recurrence relation of the form 

Pn+l(x) = X 2 p" n {x) + X(l - 2x)p' n {x) - Xp n (x), U = 0, 1, 2, . . . . 

In particular, we derive 

Po(x) = 1, Pi(x) = —x, p2(x) = 3x 2 — x, pz(x) = —15a; 3 + 15x 2 — x. 

The generating function of these polynomials is given by the series 

Pn(x) 
(2n)\ 



oo 



-2xsinh 2 (t/2) _ Pn(x) ,2n 



n=0 

Letting x = in the latter equation, we find 

p„(0) = 0, n= 1,2,.... 

The leading coefficient a n>n of these polynomials can be calculated by the formula 

a n>n = (-l) n (2n - 1)!! = (-l) n l • 3 • 5 . . . • (2n - 1), n£ff. (2.4) 

The following lemma proves an upper bound for the system p n (x). Indeed, we have 
Lemma 1. Let x > 0, n eN, e e (l - 1 



Jl + A(l-e) 2 
a e 0, 2arccos — 



and 

(2.5) 



Then 



\ P ( X )\<J (2 4 "-l)(4n)!sin(q/2) e" 
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Proof. In fact, taking representation (2.2), we apply the Schwarz inequality together with 
integral formula for Bernoulli numbers B± n , n = 1, 2, . . . [1] 

oo 4n-l /o4n _ i\ 4n 

L_ dT = -3*3-^, (2 .7) 
and relation (2.16.51.8) in [8], Vol. II. Then we deduce 

<-) i>oo r\ / ;>oo \ 1/2 

|Pn(^)| < -e^y \K iT (x)\r 2n dr < -e x I J r sinh (or) if? (a;) rfrj 

x f /°° V2 = e« f I) 2n / (l-2 4 ")sin(q/2j^: 

\J sinh(ar) / \aJ V 27m 

x (^i(2a;cos(a/2))) 1/2 , i>0,«G (0,tt). 
But in the meantime via (1.10) it is not difficult to find that 

(poo \ 1/2 

x J e- 2xcos{a/2)coshu cos\iudu 

1/2 / r oo \ 1/2 



x J e~ 2xcoB{a/2)V ^dvj < ^-z ^ y e-^rfwj = e £X ~ x 
when £ G |l - 1 and a satisfies condition (2.5). Indeed, in this case 



2 cos(a/2)vV + 1 > v + 2(1 - e) 



for any t> > since 2cos(a/2) > a/T+ 4(1 — e) 2 . Moreover, employing a sharp upper 
bound for the Bernoulli numbers — B in (see in [2]) 

d < (4n)! 



27r 4n( 2 4n-2 + 6 / 7r 2_ 1 )' 

we combine with (2.7) to arrive at inequality (2.6) and complete the proof of Lemma 
1. □ 

Remark 1. As we observe form the Stirling asymptotic formula for factorials [1], 
inequality (2.6) guarantees the absolute and uniform convergence of series (2.4) for the 
generating function on any compact set of x G [xo,A ] C M+ and t from the interval 
\t\ <t < 2/a. 
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Theorem 1. The system p n (x),x > 0, n G N can be expressed in terms of the 
absolutely convergent series 



Pn (x) = 2e x ^(-l) m m 2 "/ m (x), (2.8) 



m=l 



where I v {z) is the modified Bessel function of the third kind [5]. Moreover, we have an 
explicit formula for these polynomials 

n 

Pn(x) = J2 a k,n^ k , (2.9) 



k=l 



where the coefficients a^ n are given by 



r=0 v 7 j=0 



Finally, the following combinatorial identities hold 

k / -iX-r* / 7 \ ^ ^* 



r=0 ^ ' j=0 \ J / 



r=0 ^ 7 j=0 \ J / 

Proof. Taking the integral representation (2.2) of p n {x) and employing the formula (see 
[6]) 

K ir(x) = . . [I-ir{x] - I iT {x)] , 

2 sin(7rirj 

where I„(z) is the modified Bessel function of the third kind with its series representation 

where T(z) is Euler's gamma-function [1], we get the chain of equalities 

e x f°° (it) 211 f°° (ir) 2n 

Pn(x) = tt / • 7 • x [-f-w-fa) - I iT (x)] dr = e x / . , I- ir (x) dr 
2 J.^ sm(mr) J_ oc sm(Trzr) 



x moo . . 

- / ^-^(l + *)r(i -s)J2 

1 J-i<x> 



, *!r(* + 1 - S ) 
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The change of the order of integration and summation in the right-hand side of the latter 
equality is indeed possible by Fubini's theorem owing to the absolute convergence for each 
x > 0, namely 



Jin- 



k=0 



(:r/2) 



2k-s 



k\Y(k + 1 - s) 



ds 



< oo. 



Therefore, 



Jin- 



. 1 r(i + s)r(i- s ) fx 



r(k + i- s) 



(f) 



ds. 



(2.14) 



But the integral in (2.14) can be calculated via the Slater theorem (see, for instance, in 
[8], Vol. Ill) after 2n — 1 times differentiation with respect to z under the integral sign. 
This operation is allowed by virtue of the absolute and uniform convergence. Thus we 
obtain 



/' 

j — 1 



Jin- 



^r(i + s )r(i- s ) 



r(k + i-s) 



d \ 



(2n-l) -ioo 



r(i + 3 )r(i- 3 ; 
r(fc + i - s) 



2 " 1 2 s) 



(2n-l) oo 



E 



;-i) m+1 z m+1 r(2 + m) 



" Q ml r(2 + fc + m) 



2 " 1 *£) 



(2n-l) oo 



E 



-1)"V" m 
^ (m + /c)! 



Hence differentiating under the series sign due to the absolute and uniform convergence, 
letting z = x/2 and combining with (2.13), (2.14), we derive the representation 



fc=0 



m=l 



m=l 



where the inversion of the summation order is motivated by the absolute convergence of 
the iterated series. Thus we proved (2.8). 

On the other hand, calling relation (5.8.5.3) in [8], Vol. II, it is not difficult to find via 
differentiation with respect to a parameter owing to the absolute and uniform convergence 
that 

2 V m 2n I m (x) = lim — — e xcosha . 
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We calculate the 2n-th derivative in the right-hand side of the latter equality employing 
the Hoppe formula [5]. Precisely it gives, 



cosh a = e xcoshaSp x k ( COslia) fc 3 (f n ^ 

da 2n ^ ^ j\(k-j)\ da 2n 

= .cosh a V T fc V (- COshQ ) fc " J V (A e^~l> 

2^ 2ij\(k-j)\ da? n ^\r 

k=0 j=0 J v J ' r=0 v 7 

= „•«">«£ * f: ( ~™^ E (i) (2r)—»(2 r - i)»eP-*. 

fc=0 j=0 J '' r=0 ^ 7 

Therefore, passing to the limit when a — > 0, we find the identity 

2 f m 2 "U*) == «■ t E (*) ^ E (j) (* - ^ 

m=l fc=0 ' i=0 ^ 7 r=0 ^ 7 



(-i)* ^ (-1)' fk\ ^ (-1)' ft-rv, .. 

e E— E— (J E— ( j jo-* 

fc=0 r=0 v 7 j=0 V J 



2n 



On the other hand, since (— l) m I m (x) = I m (—x), we return to (2.15) and appealing to 
the right-hand side of the latter equality, we derive 



In , k / i w /i \ k—r 



^e^w: E^er>-^" 

fc=0 r=0 v 7 j=0 

However, as it is proved in [11], p„(x) is a polynomial of degree n. Therefore, all coefficients 
in front of powers x k , k = n + 1, n + 2, 2n are surely equal to zero. Thus we 
establish the explicit formula (2.9) with coefficients (2.10) and combinatorial equality 

(2.11) . Taking into account the value (2.4) of the leading coefficient a„ jn , we get identity 

(2.12) and complete the proof of Theorem 1. □ 

Finally, in this section we will established an exact asymptotic behavior of p n (x), when 
n — > oo and x > is a fixed number. Precisely, we have 

Theorem 2. Let x > and j3 G (0, 7r/2) be fixed numbers. Then 

/x 6x(-l) n sin/3 (2n)\ T / „ f 1\\ 
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Proof. Indeed, employing representation (2.2) and choosing a fixed parameter j3 G (0, 7r/2), 
we write 

Pn ( x ) = A >- e x / T 2n e-P T [cosh(pY) + sinh(/3r)] K lT (x) dr. 

n Jo 

Hence the Parseval identities for the cosine and sine Fourier transforms [9] together with 
relations (2.5.31.4) in [8], Vol. I and (2.16.48.20) in [8], Vol. II drive us to the equalities 

Pn(x) = -e x / — — — , [cos ((2ra + l)arctg(y//3)) cos(a; sm(/3) sinh y) 

Jo (P + y 2 ) n+L ' 2 

+ sin ((2n + l)arctg(y//?)) sin(x sm(/3) sinhy)] dy 

2( — l) n (2n)' f°° g-zcos(/3)cosh?/ 

= : e x ^ 2 ~ y2 )n+i/2 cos ^ 2n + 1 ) arct g(f// 3 ) - xsm(/3) sinhy] cfy 

-(2n + l) (log(v^T7) 

+iaxctg(y/P))]dy]. (2.17) 

The asymptotic behavior for large n of the latter integral in brackets can be treated by 
the Laplace method [7], Ch. 4. Indeed, it has the form 



Re 



e -xcos(p+iy) exp 



POO 

/ q(y)e-^ +1 ^dy, 
Jo 



where p(y) = log(/3 + iy), q(y) = e xcos (P+ i y). Thus the integral 

where the calculation of three first coefficients a s will be enough to determine the main 
term of the expansion of its real part. Namely, using formulas for coefficients in [7], we 
have 

a = -ip e - xcos P, ai = -i(3(l+x(3sm(3) e~ xcos ^ 
a 2 = -x(3 sin f3 - l A e -*™ s p ^ x ^ cos ^ + x gin 2 + 
Therefore, returning to (2.17), we find 

. 6x(-lfsin/3 (2n)\ / 
^ = ^(2n + l)3, 6 t, 1 + ° UJ J ' » - °°' 

which proves (2.17). 

□ 
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3 Representations of the Yor integral 

This section will complete our goal to represent integral (1.1) in a different form. Namely, 
we will start relating F t (r) with the inverse Kontorovich-Lebedev transform (1.4). Indeed, 
taking (1.4) with v = it and integrating by parts we come out with the representation of 
the modified Bessel function 



.(r) = r f 
Jo 



-rcoshu 



sinh-u sin rudu, r > 0. 



(3.1) 



Hence applying the Parseval equality for the sine Fourier transform [9] with the relation 
(2.5.36.1) in [8], Vol. 1, integral (1.1) becomes (r,t > 0) 



F t (r) = 



TT 2 /2t rOO 



f 

Jo 



rK lT {r) 



exp 



t(r - n/tf 



cxp 



t(r + ir/ty 



dr 



2 f°° _t 2 
= — - / e _ 2 r r sinh 7iTK iT (r)dT. 
rir 2 J 

Hence form (1.2) we obtain reciprocally for each t > 

POO 

/ K tT (r)F t (r)dr = e-^ 2 
Jo 

and via Parseval equality (1.3) we derive the value of the integral 

poo 2 r°° 

/ \F t (r)\ 2 rdr = — rsinh(7rr) e~ tT * dr 
Jo 7T J 



(3.2) 



(3.3) 



2e w2/4t [°° 

J-OO 



Recently (see [12]), the author introduced the following heat kernel for the Kontorovich- 
Lebedev transform 



h t (x,y) 



2 r°° 2 

— - / e~ tT I 2 t sinh txt K iT (x)K ir (y)dT. 

X7l J 



(3.4) 



Hence appealing to the Macdonald formula (1.15) and Fubini's theorem to interchange 
the order of integration, we come out with the representation of (3.4) as a translation 
operator of the Yor integral for convolution (1.12). Precisely, minding (3.2) it gives 



h t (x, v) = — 



F t (r)dr, x,y > 0. 



(3.5) 
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Meanwhile differential and convolution properties of the Yor integral (1.1) are given 
by the following 

Theorem 3. The function F t (r) is infinitely differentiable of variables (r,t) G R+ x 
R +; satisfying the estimate 



d m F t {r) 



nl/4— m p tt 2 /t 

< tm+ 3/ 47r n/s K^(2r)T^(Am + 5/2), m G N„. 



(3.6) 



Moreover, F t (r) is a solution of the generalized diffusion equation (u = u(t,r)) 



du 9 d 2 u du 
dt dr 2 



_l_ r r u 

or 



and satisfies the index law in terms of convolution (1.12) 

F t (r) = (F tl *F t2 )(r), h + t 2 = t. 
Finally, for the Yor integral the following integral equation takes place 

ht 1 (r,y)F t2 {y)dy = F tl+ta (r). 



I 

Jo 



(3.7) 



(3.f 



(3.9) 



Proof In fact, employing the following inequality (see [11], [13]) for derivatives of the 
modified Bessel function with respect to x 



d m KiJr) 



dx r 



< e- 5T K m (r cos S), x > 0, r > 0, 5 E 0; |) 



0;-), m = 0,l, 



(3.10) 



it is not difficult to verify that all positive t the integral in the right-hand side of the 
latter equality in (3.2) together with derivatives of any order with respect to r converge 
absolutely and uniformly by r > ro > 0. Therefore F t (r) is infinitely differentiable with 
respect to r. Similar motivation can be done for the derivatives of the order m G No with 
respect to t > and it gives the expression 

d m F t (r) 2 1 -" t (-l) 



dt r 



i/ i \m roc 

e -^ 2 /2 r 2m+i sirih7rr K iT {r)dT, m = 0,l, ... . (3.11) 
n Jo 



Hence the Schwarz inequality and relation (2.16.52.6) in [8], Vol. II yield 



d m F t (r) 



)1— m 



< 



< 



dt m 

nl/2-m 
2 ^V2 



71 z 



K 1 (2r) 



e -tr 2 +27vr T 2(2m+l) dT 
1/4 



1/2 



\ 1/2 

Kl(r)drj 



-tr 



+AnT dl 



e -tr» r 4(2m+l) dr \ 
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2l/4-r. 



e 7r2/i A' 1/2 (2r)r 1 / 4 (4m + 5/2). 



^-m+3/4^11/8 

Thus we proved (3.6). 

Further, as it follows from (1.5), (1.6) and absolute and uniform convergence of the 
corresponding integrals, formula (3.11) can be written as the following partial differential 
equation 

d m FAr) 

where m — 0, 1, . . . and A™ is m-th iterates of the operator (1.5). In particular, letting 
m = 1 we obtain that the Yor integral (1.1) satisfies the generalized diffusion equation 
(3.7). 

Next, equality (3.8) is a direct consequence of (3.3), factorization equality (1.14) and 
the uniqueness property of the Kontorovich-Lebedev transform [11]. In order to prove 
(3.9), we multiply h tl (x,y) by F t2 (y) and integrate by y over R + . Then using (3.3), 
(3.4) and the Fubini theorem, which is applicable due to the absolute convergence of the 
iterated integral (it can be verified, employing inequalities (3.6) with m = and (3.10)), 
we get the result. □ 

Finally, we will establish a representation of the Yor integral (1.1) in terms of polyno- 
mials (2.1) and the heat kernel (3.4). 
We have 

Theorem 4. Let r, t > 0. Then the Yor integral satisfies the following equation 
where * denotes convolution (1.12). Moreover, it can be rewritten in the form 



where 



f°° _ du 
a k (r,t)= / e~ u ht(r,u)p k (u) — 

Jo 



u 



and h t is the heat kernel (3.4). 



Proof. Calling again integral (3.2), we expand the hyperbolic sine in Taylor's series. Then 
changing the order of integration and summation due to the estimate 

poo 00 /__\2fc-l 00 ~2fc-l /-oo 
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Mr) ^ (*y/yi) r(fc + i/2) 

1^ Ok-lV <00, 



vrv^tt (2k-l)\ 
we come out with the representation 

W = ^E(^/ e"-/ 2 T ^, r(r)dT . (3 , 6) 

But the integral in (3.16) can be treated with the convolution for the Kontorovich-Lebedev 
transform (1.12), its factorization property (1.14) and inversion formula (1.4). Thus 
minding formulas (2.2), (3.2) and (3.3), we derive 

— e-"'' 2 T 2k K iT (r)dT = ( e - r ^ * F t (r)) . 

rir 2 Jo it \ r J 

Substituting in (3.16), it gives 

and equality (3.14) is proved. On the other hand, employing (3.5), one can express the 
latter convolution in terms of the heat kernel (3.4). Precisely, we find for each fceN (see 

(1-12)) 

which proves (3.15). 

□ 
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